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ABSTRACT 


Jlo  study  granular  jumps  generated  in  the  gravity  flow  of  granular 
materials  down  an  inclined  plane.  The  equations  governing  the  granular  flow 
are  reduced  to  a  sequence  of  boundary  value  problems  of  linear  ordinary 
differential  equations  by  means  of  an  asymptotic  method.  Their  solutions  are 
used  to  determine  the  Burgers  equation ,  which  possesses  a  progressive  wave 
solution  to  describe  a  smooth  granular  jump.  The  wave  speed  and  a  criterion 
for  the  stability  of  the  granular  jump  are  also  obtained. 
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SIGNIFICANCE  AND  EXPLANATION 


Problems  of  flow  of  granular  materials  have  attracted  much  attention 
recently  because  of  their  importance  in  industrial  applications  and 
geophysical  situations.  It  is  observed  in  laboratory  tests  that  a  granular 
jump  may  appear  in  the  gravity  flow  of  a  granular  material  down  an  inclined 
plane.  The  granular  jump  may  be  considered  as  a  wave  front  moving  with 
constant  velocity  on  the  free  surface  of  the  granular  material  and  connecting 
one  uniform  depth  to  another.  The  main  contribution  of  this  report  is  to 
develop  an  asymptotic  method,  which  can  describe  the  speed,  wave  profile  and 
stability  of  the  granular  jump. 
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GRANULAR  JUMPS  DOWN  AN  INCLINED  PLANE 


M.  C.  Shen 


1.  Introduction 

In  recent  years  there  has  been  much  interest  in  the  study  of  the  gravity 
flow  of  granular  materials.  Needless  to  say,  granular  flow  problems  occur 
frequently  in  industry  as  well  as  in  nature  such  as  handling  of  grains  and 
sands,  snow  avalanches  and  rock  falls.  It  is  reported  in  laboratory  tests 
(Savage,  1979)  that  smooth  granular  jumps  may  be  observed  when  granular 
materials  are  moving  down  an  inclined  chute.  As  seen  from  loose  snow 
avalanche  or  earth  slide  in  nature,  such  a  granular  jump  also  seems  to 
appear.  On  the  other  hand,  since  Goodman  and  Cowin  (1972)  proposed  their 
continuum  theory  for  flow  of  granular  materials,  only  a  few  analytical  results 
based  upon  their  equations  have  been  found,  for  example,  steady  state 
solutions  by  Goodman  and  Cowin  (1971),  and  solutions  of  the  linearized 
equations  by  Nunziato  and  Walsh  (1977).  The  main  purpose  of  this  report, 
therefore,  is  to  develop  an  asymptotic  method  for  the  study  of  the  time- 
dependent  granular  flow  under  gravity  down  an  inclined  plane,  which  provides  a 
means  to  describe  the  development  of  a  smooth  granular  jump.  Some  results 
concerning  granular  jumps  have  been  obtained  by  Morrison  and  Richmond  (1976), 
and  Savage  (1979).  However,  their  approach  is  based  upon  algebraic 
considerations  of  conservation  laws  for  hydraulic  jumps. 
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In  this  report  we  adopt  the  equations  formulated  by  Goodman  and  Cowin 
(1972)  for  the  so-called  cohesionless  Couloetb  granular  material,  and  our 
approach  used  here  is  essentially  an  extension  of  the  method  we  developed  in  a 
previous  report  (Shen,  1983)  for  the  compressible  viscous  flow  down  an 
inclined  channel.  We  briefly  explain  the  method  as  follows.  Assume  that  the 
granular  flow  down  an  inclined  plane  is  two  dimensional.  We  choose  a 
coordinate  system  moving  with  the  wave  front  of  a  granular  jump  so  that  the 
horizontal  axis  is  parallel  to  the  plane.  A  small  parameter  t  as  the  ratio 
of  the  vertical  length  scale  H  to  the  horizontal  length  scale  L  is 
introduced.  We  assume  that  the  vertical  velocity  component  is  much  smaller 
than  the  horizontal  velocity  component,  and  consider  the  granular  flow  for 
large  time.  Then  we  use  e  to  stretch  various  variables  non-dimenionalized 
by  appropriate  units  so  that  these  assumptions  are  met.  The  solution  of  the 
governing  equations  is  expanded  in  an  asymptotic  series  in  integral  powers  of 
c,  and  the  equations  for  the  successive  approximations  determine  the  wave 
speed  and  an  evolution  equation  for  the  free  surface,  which  is  the  well-known 
Burgers  equation.  The  signs  of  some  coefficients  in  the  Burgers  equation 
determine  whether  the  equation  is  ill-posed  and  provide  a  criterion  for  the 
stability  of  the  flow. 

The  governing  equations  used  here  involve  four  material  constants;  two  of 
them  are  related  to  the  equilibrium  stress  and  the  other  two  appear  in  the 
dynamic  stress.  We  first  develop  a  general  method  in  which  all  the  material 
constants  are  assumed  to  be  finite.  However,  the  calculations  of  the  wave 
speed  and  the  coefficients  in  the  Burgers  equations  are  rather  prohibitive. 

To  simplify  the  derivations,  we  shall  consider  a  special  case,  that  is,  we 
choose  one  of  the  material  constants  to  be  large  in  comparison  with  the 
others.  In  this  case  the  equilibrium  solution  reduces  to  the  one  for  the 
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incompressible  viscous  flow  down  an  inclined  plane,  and  the  wave  speed  and  the 
coefficients  in  the  Burgers  equation  can  be  explicitly  expressed  in  tens  of 
the  regaining  constants.  He  also  note  that  many  other  constitutive  equations 
have  been  proposed  for  granular  materials  lately  (Cowin  and  Satake,  1978 j 
Savage,  1979|  Nunziato,  Passman  and  Thomas,  1980;  Savage  and  Jeffrey,  1981; 
Sayed  and  Savage,  1983).  Most  of  them  are  generalizations  of  the  original 
formulation  due  to  Goodman  and  Cowin  (1972).  The  basic  ideas  developed  here 
could  also  apply  to  these  more  general  cases. 

We  formulate  the  problem  for  the  general  case  in  Section  2.  The  wave 
speed  and  the  Burgers  equation  are  derived  in  Section  3.  In  Section  4,  we 
consider  the  special  case.  In  Section  5,  we  define  the  criterion  of 
stability,  derive  an  expression  for  the  wave  profile  of  a  granular  jump  and 
present  some  problems  for  further  study. 


2.  Formulation 

He  consider  a  cohesionless  Coulomb  granular  naterial  flowing  down  a  rough 
inclined  plane  so  that  we  My  impose  the  no-slip  condition  there.  The  flow  is 
assumed  to  be  two  dimensional  and  a  coordinate  system  moving  at  a  speed  o* 
in  the  x*-direction  is  chosen  so  that  the  x*-axis  is  parallel  to  the  plane  and 


the  y*-axis  is  positive  upward  (Figure  1).  The  governing  equations  are  the 


following: 


V*#  ♦  7*  •  (V*q* )  -  0  , 

yv*(q*#  +  q*  •  7*q*)  ■  7*  •  T*  +  yv#q* 


subject  to  the  boundary  conditions: 


At  the  free  surface  £< 


at  the  rough  rigid  plane. 


-  n*(x*,t*)  -  o  , 

T*  •  n*  -  0  , 


h*  •  n*  =  0  , 


n*  ♦  u*n*  -  v*  *  0  ; 

t*  x* 


q*  _  -o*  v*  »  v*  . 
H  P 


Here  7*  »  (3/3x*,3/3y* )  *  (3/3x*,3/3x*) ,  v*  is  the  volume  distribution 
function,  Y  is  the  mass  density  assumed  to  be  constant,  q*  -  (u*,v*)  - 
(u*,u*>  is  the  velocity  of  the  flow,  T*  is  the  stress  tensor  defined  by 

T*  -  I-B*(v*)2  +  a*(7*v*)2  +  2o*v*(7*)2v*] 5 .  -  2a*v*  v*  + 

ij  x*  x* 

X*7*  •  q*6 ,  +  211*0?,  , 

JLj  Aj 

o*  *  (u*  *  +  u*  *)/2  , 

■ij  jx* 

g  ■  (g  sin  0,  -g  cos  9)  is  the  constant  gravitational  acceleration,  9  is 
the  angle  of  inclination  of  the  plane  to  the  horizontal  assumed  to  be  greater 
than  the  angle  of  internal  friction  but  less  than  ir/2,  h*  is  the  equi¬ 
librated  stress  vector  given  by  h*  “  2a*7*v*,  and  n*  is  the  vector 


(n*#,-1)  normal  to  the  free  surface*  For  simplicity,  we  assume  that  a*, 
0*,  u*  and  v*  are  positive  constants,  and  X*  is  also  constant  but 
satisfies  3X*  +  2v*  >  0. 

We  introduce  the  following  nondimensional  variables: 


v  -  v*/v*,  (u, v)  -  (u*,£v*MgH)  , 

P 

t  -  e"2t*(H/g)~1/2,  (x,y)  *  (E'V.yMH-1  , 

a  -  o*(gH)_1,  $  =  0*(yv*gH)  1  , 

P 

a  -  aMyAgH3)'1,  <X,v)  *  <X*,y*)  (yv*(gH)  1/2H)~1 

P 

e  «  H/L  . 


In  terms  of  them,  (1)  to  (6)  become 

ev  +  (vu)  +  (vv)  =  0  , 

t  x  y 

2  2 
ev(e  u  +  uu  +  vu  )  *  -20cw  +  2ave(e  v  +  v  ) 

t  x  y  x  xx  yy  x 

2  2 

+  (X  +  u)e  (u  +  v  )  +  u(e  u  +  u  )  +  v  sin  9  , 

x  y  x  xx  yy 

2  2 
e  v(ev  +  uv  +  w  )  =  -20vv  +  2av(c  v  +  v  ) 

t  x  y  y  xx  yy  y 

2 

+  (X  +  jj)e(u  +  v  )  +  ye(e  v  +  v  )  -  v  cos  0  ; 

x  y  y  xx  yy 

at  y  *  n(x,t) , 

2  2  2  2  2  2  2 
t-Bv  +  a(e  v  +  v  }  +  2av(£  v  +  v  )  -  2ae  (v  ) 
x  y  xx  yy  x 

2 

♦  Xe(u  +  v  )  ♦  2wcu  ]cn  +  2aev  v  -p(u  +ev)=0  , 

xy  xx  x  y  y  x 

2  2  2  2  2 
[-2aev  v  ♦  g(u  *  e  v  )]En  +  Bv  -  a(E  v  +  v  ) 
x  y  y  x  x  x  y 

2  2 

-2av(£  v  +  v  )  ♦  2o(v  )  -  Xe(u  +  v  )  -  2pev  =  0  , 

xx  yy  y  x  y  y 

e2t)  v  -  v  =0  , 

xx  y 
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3.  Wave  speed  and  the  Burgers  equation 

Assume  that  u,  v,  v  and  n  possess  an  asymptotic  expansion  of  the  form 


♦  “  +  e4i  +  e  V" 


(15) 


Without  loss  of  generality,  we  may  assume  a  *  oQ  +  eo^.  Substitution  of  (15) 
in  (7)  to  (14)  yields  a  sequence  of  equations  and  boundary  conditions  for  the 
successive  approximations.  The  equations  for  the  zeroth  approximation  are 


at  y  =  0, 


at  y  =  -1, 


V  u 


Oyy 


-Vq  sin  0 


2  a  v„  -  =  cos  0  ; 

Oyyy  Oy 


u0y 


=  0 


2  a  v„  -  0v„  =  0 
Oyy  0 

V  =  0  i 
oy 


-a.,  =  1 

0  0 


(16) 

(17) 

(18) 

(19) 

(20) 

(21) 


where  u_  is  assumed  to  be  a  function  of  y  only  and  v 


0.  The 


solutions  of  (16)  to  (21)  for  uQ  and  can  be  easily  found,  and  a 

discussion  of  these  solutions  is  given  in  Goodman  and  Cowin  (1971).  Following 
them,  we  define 


l  =  (P/a) 


1/2 


and  M  =  cos  0/(20) 


(22) 


and  obtain 


Vg  =  Ml  sinh  £y 


+  (1  -  M  +  Mi  ^  sinh  t)(1  +  cosh  %)  ^(1  +  cosh  Jly)  -  My 


(23) 


U0  *  *0  “  °0 

=  u  1  sin  0{-Ml  3(sinh  l  +  sinh  iy)  +  (1-M+M£  1  sinh  i)x 

(1  +  cosh  t)  1 1 ( 1  — y2 )/2  +  l  2(cosh  %  -  cosh  ty)]  +  (M/6)(1+y3) 

+  Ml  2 ( 1+y ) }  -  Op. 


(24) 
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For  small  l,  we  may  expand  vQ  and  in  integral  powers  of  l  and 

obtain 


V0;  1  +  t2[(M/6)y3  +  y2/4  +  (M/6  -  1/4)]  . .  , 

<!>0  2  M*1  sin  0{  ( 1~y2  )/2  +  f2  [  ( 1  ~y2  )  (M/12  -  1/8)  +  (1-y4)/48 

-  M(1+y5)/120]+ . }  . 


(25) 


(26) 


We  see  from  (25)  and  (26)  that  as  l  *  0,  vQ,  uQ  tend  to  the  solutions  for  an 
incompressible  viscous  flow  fown  an  inclined  plane  (Landau  and  Lifshitz, 

1959).  This  result  motivates  our  choice  of  a  as  a  large  parameter  in  the 
next  section  so  that  expressions  for  the  wave  speed  and  the  coefficients  in 
the  Burgers  equation  can  be  greatly  simplified.  On  the  other  hand,  if 
l  *  ®,  vfl  tends  to  a  discontinuous  solution  when  M  /  1  and  becomes  zero 
at  y  =  0  (Goodman  and  Cowin,  1971).  Therefore,  we  shall  not  study  this  case 
in  the  sequel. 

\ 

The  equations  for  the  first  approximation  are 

0  , 


at  y  *  0, 


a(v  v 


(v  u.  +  vnu.)  +  (v.v,) 
10  Olx  Oly 

Vi  u,  =  -v ,  sin  9 

iyy  i 

a  v.  -  8  v,  =0 
lyyy  1y 


u  +  u  n  =  o 

1y  Oyy  1 
+  __  v.)  -  8  v„v,  =  -a  v„v 


0  Iyy  Oyy  1 


0  1 


0  Oyyy  1  ' 


at  y  =  -1, 


v,  +  v  n  «  o 

1y  Oyy  1 

u  n  “v  *  0 
u0  lx  1 


“V  vi =  vi 


As  observed  from  (28)  to  (32)  and  (34),  we  may  express  vt  and  u1  as 


V1  =  Vi'  U1  "  Vi  -  °1 


(27) 

(28) 

(29) 

(30) 
(31  ) 

(32) 

(33) 

(34) 

(35) 
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at  y  -  0 


2  2 
<-»V2*Voy,,nix  '  ""V“.yVu0„,’2  *  “Oyyy"/21  ‘  0 

»'2WnV2  *  V1  *  2V(l“lyni  *  2Vj>  *  2,(VOyy,2"l  *  % 

2 

+  2\>  v  nJ  -  2<xIva(va  **  ♦  vo  n-/2  +  v  n  +  v_ 
Oyy  iy  i  o  Oyyy  2  Oyyyy  i  iyyy  1  2; 

♦  VA  (VA  o^/2  *  v4  o,  ♦  vo  ♦  V4V«  )  ♦  vjv  n,  ♦  V, 
Oyy  Oyy  1  1y  1  2  1  lyy  1  Oyyy  1  Iy 


Mu.  ♦  v.  )  -2u  v-  *  0  , 

lx  iy  1y 


V2y  +  V1yn1  +  V0yyV2  =  °  ' 

nit  +  V2x  +  Vlx  *  UOyVlx  "  V2  "  V1yni  “  °  ' 


at  y  =»  -1, 


u.*  =  v_  «  v_  *  0 
2  2  2 


From  (17),  (19),  (24),  (35),  (36)  and  (42),  we  may  rewrite  ( 
and  (53)  as  follows: 

S<"2’  ■  “  »2yy  "  "U2  ,in  9  *  'VVVS  *  2BVl  ‘  2aV 
*  VlV”„  ' 

W  ■  2aV0V2yyy  '  a6V0V2y  '  * 1  IXt“>  (*.y*f  1yy>  * 
at  y  =  0, 

W  “  UU2y  =  (_RV0+2aV0V0yy)T11x  "  U  (*  1y+Vyy/2  ^  "  U*0yy 


B2(v2)  *  2av0v2yy  -  Pv0v2  =  f*(*1+f1y>  -  2Wf1y)»)1x 

+  {6<V0yy+*1*2Vly)  "  2af(V0yy)2/2  *  <*1y,2/2  *  V' 


♦  "W  *  VyVlyy  +  *1(V0yyy  +  *1yy)nni 


2av„v„  n.  . 

0  Oyyy  2 
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(58) 
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Oyy 


In  principle,  we  can  solve  (62)  to  (67)  for  and  and  in  turn  u2 

and  can  be  determined  from  (60),  (61).  Now  we  are  in  a  position  to 

determine  the  Burgers  equation.  We  integrate  (46)  with  respect  to  y  from 
y  *  -1  to  y  *  0  and  make  use  of  (18),  (33)  to  (35),  (52)  to  (53),  (60)  and 
(61),  to  obtain 


/!i  <*  v  +  /V<Vix  +  *3ni  +  W(W  +  Wr°i)T,i 


+  V0(*2nix  +  *3ni  +  *4n2J1xdy  +  ^1 (0  J  ^0 (°  >“°01 ni ni: 


(68) 


+  nit  +  (Wn2x  +  (*lnr°1)nix  +  ‘f1y(0,ninix  “  0  * 


Note  that  the  coefficient  of  Hj  vanishes  because  of  (45).  By  rearranging 
the  terms,  we  finally  obtain  the  Burgers  equation 

Vlt+"ln1x+  Vln1x"  Vlxx  '  {69) 

Which  is  our  main  result.  Here 

m0  =  /-I  *1dy  +  1  '  m1  *  “°im0  ' 

m2  -  2  /®t  N^VV  +  +,♦,  +  V0*3)dy  +  *1(0>[V0>-V  +  V0)  "  f1y(0)  ' 

m3  =  “  /VVW  +  V2]dy  * 


(69)  can  be  used  to  study  the  development  of  a  granular  jump. 
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4.  A  special  case 

In  Section  3,  we  derived  the  expressions  for  the  wave  speed  oQ  and  the 
coefficients  for  the  Burgers  equation.  However,  they  depend  upon  6,  a,  0,  X 
and  w,  and  are  too  complicated  for  applications.  Therefore,  as  Motivated  by 
the  discussion  of  vQ  and  given  in  the  last  Section,  we  assume  a  *  1/e 

to  be  a  large  parameter.  Substitution  of  (15)  in  (8)  to  (14)  also  yields  a 


sequence  of  equations  and  bounding  conditions.  Hie  equations  for  the  zeroth 
approximation  are  the  following: 


Mu.  =  -v  sin  0  , 

Oyy  0 

(70) 

V  »  0  ; 

Oyyy 

(71  ) 

at  y  -  0, 

u0y  *  0  *  (72) 

v  *  0  ,  V-  =  0  ; 

Oyy  Oy 

(73) 

at  y  =  -1, 

u„  *  -O  „ ,  V  =  1 

0  0  0 

(74) 

It  is  easily  seen 

that 

vo  *  1  ' 

(75) 

2 

uQ  *  (sin  0/(2m)J(1  ~  y  )  ~  ao 

(76) 

(75)  and  (76)  are 

solutions  for  the  incompressible  viscous  flow  down  an 

inclined  plane. 

The  equations 

for 

the  first  approximation  are 

,viuo  *  Vi’x  *  ‘ViV  "  0  ' 

(77) 

M  u,  *  -2v,  -  v,  sin  0  ; 

lyy  lyyx  1 

(78) 

2v  -  cos  0  , 

lyyy 

(79) 

at  y  ■  0, 

U1y  +  U0yyni  =  °  ' 

(80) 

v  =  0  , 

lyy 

(81  ) 
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0 


(82) 


viy  “ 

u.n,  “  v, 
0  lx  1 


0  j 


(83) 


at  y 


■1. 


U1  = 


V  V1 


(84) 


(85) 


It  is  found  from  (78)  to  (82),  and  (84)  that 

*  (cos  0/12 ) ( 1  ♦  y3) 

2  5 

u  |  ■  [sin  0  cos  0/(2^v )]  (9/10  -  y  -  y  /10)  -  ♦  (sin  B/|i )  ( 1*y  )n  ^  .  (86) 

We  now  integrate  (77)  with  respect  to  y  from  y  =  -1  to  y  and  obtain,  by 
(75),  (76),  (85)  and  (86), 

v1  *  u1x<Sy  *  -tain  0/(2w))  (1+y)2t)1x  .  (87) 

It  follows  from  (76),  (83)  and  (87),  that 

[sin  0/(2p)  -  <J0]h1x  -  -sin  0/(2|i)n1x  . 


Assume  that  n.  Z  0  and  we  have 
lx 

aQ  *  sin  0/y  . 

Now  we  proceed  to  the  equations  for  the  second  approximation 

vit  *  ‘Vo  *  Vi  ♦  voV,  *  IV1V1  *  w,  ■  0  • 

Vi,  *  Voy  ”  *  “°2yy  *  *2  *1“  9  1 

-  20v  +  2 (V  +  v  v  )  -  v  cos  0=0  » 

1y  2yyy  1  lyyy  1 


at  y  -  0, 


<-0+2v,  )n,  -  u(u.  +u,  n,+u.  n.+u.  n,/2)  =  o 

lyy  lx  2y  1y  1  Oyy  2  Oyyy  1 

2 

20v  =  )v  )  -  2(v  -v,v  +v,  n,)  =  0  , 

1  1y  2yy  1  lyy  lyyy  1 

v.  +  v.  n,  *  o  , 

2y  lyy  1 


V  +  V2x  *  Vlx  +  U0yVlx 


V2  -  V1yni 


(88) 

(89) 

(90) 
(91  ) 

(92) 

(93) 

(94) 

(95) 


at  y  -  -1, 


u_  -  v_  =  v_  =  0 
2  2  2 

From  (85),  (91),  (93),  (94)  and  (96),  we  find  that 

2  5  2 

■  (0  cos  0/24) (-9/10  ♦  y  +  y  /10)  +  (cos  9/4)(1-y  )n^ 


(96) 


(97) 
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Then  from  (76),  (86),  (87),  (90),  (92),  (96)  and  (97)  it  is  obtained  that 
yu2  ■  (-6  -  cos  8/2  -  sin28/(8y2)  -  By  +  y2cos  8/2  +  sin28(-y2/6  +  y*/12)]n1x 

+  cos  8  sin  8(5-6y2+y4)n^/48  -  sin  8(1+y)n2  +  sin  8(1+y)ri2  (98) 

-  B  cos  8  sin  8(155-189y2  +  35y4  +  y?)/10080  . 


Making  use  of  (74),  (76),  (83)  to  (88)  and  (95)  to  (98),  we  integrate  (89) 
with  respect  to  y  from  y  =■  -1  to  y  =  0  to  obtain 


It  1  lx  2  1  lx  3  Ixx 

where  the  coefficient  of  n2  vanishes  because  of  (88),  and 


m  =  cos  8  sin  8/{24p)  -  a,  , 


m2  -  sin  8/y  , 

m3  -  B/2  +  cos  8/3  +  23  sin28/(240u2)  . 


(99) 
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S.  Discussion 


It  is  trail  known  that  the  Burgers  equation  as  given  by  (69)  or  (99)  is 
ill-posed  if  the  coefficients  mg  and  m3  are  of  opposite  signs.  Therefore, 
we  may  use  this  result  to  define  the  criterion  of  stability.  Since  in  (69), 
and  depend  upon  9,  a,  6,  X  and  u,  we  say  a  granular  jump  is 

stable  if  mg  and  a3  are  of  same  sign  i.e.  m^/nig  »  f  (6  ,a,0  ,X  ,p )  >  0. 

For  (99),  nig  *  1  and  m^  >  0.  Hence  a  granular  jump  is  always  stable  in 
this  sense. 

The  solution  method  for  the  Burgers  equation  is  also  trail  known  and  tra 
only  give  a  brief  discussion  of  the  steady  state  solutions  of  (69)  and  (99)  in 
reference  to  a  moving  frame.  Assume  all  the  coefficients  in  (69)  are  not 
zero,  Bj/Bq  >  0  and  *  0  aa  x  +  "■  It  is  easily  found  that 

n1  *  -(n^/m-HI  +  tanh(m^  (x+c)/(2m3))) 

where  c  is  an  arbitrary  constant  and  tra  require  m^/mj  <  0.  Since  in  (69), 
m3  =  -BgO^,  Bj/nig  >  0  and  <  0  i«ply  >  0.  Since 

a  ■  Og  +  eo^  >  Og,  we  say  the  granular  jump  moves  at  a  supercritical  speed. 

If  Bg/m2  >  0,  then  n.j  ♦  2mgC1/m2  >0  as  x  ♦  In  this  case  n1 

decreases  smoothly  from  2m^/a to  zero  as  x  increases.  On  the  other 
hand,  if  mg/mj  <  0,  then  increases  smoothly  from  u>0a3/m2  <0  to  zero 

as  x  increases.  For  the  special  case  considered  in  Section  4,  mQ  >  1, 
m2  >  0,  m3  >  0.  Hence  m^/m^  <  0  implies  >  cos  6  sin  0/(24u)  >  0  and 
the  speed  of  the  granular  jump  is  supercritical.  Furthermore,  m2  >  0 
implies  that  we  always  have  a  decreasing  granular  jump. 

Finally  we  consider  some  problems  for  further  study.  If  the  free  surface 
of  the  granular  material  is  replaced  by  a  rigid  plane,  then  one  boundary 
condition  is  lost  and  we  do  not  have  enough  boundary  conditions  to  determine 
the  flow.  In  fact,  in  the  gravity  flow  of  a  granular  material  between  two 
vertical  plates  (Goodman  and  Cowin,  1971),  a  slug  flow  aopears  in  the  central 
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part  of  tha  vartical  channel.  It  ahould  be  of  interest  to  extend  the  approach 
to  the  study  of  the  granular  flow  between  two  parallel  inclined  planes. 
Furthermore,  if  mg/m^  <  **»  investigation  of  the  ill-posed  problem  is 

certainly  of  importance.  Finally  the  asymptotic  method  used  here  is  based 
upon  the  existence  of  a  solution  given  by  (IS).  A  justification  of  the 
method,  needless  to  say,  will  he  a  significant  contribution  to  granular  flow 
problems . 
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